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ABSTRACT 

Context. Data from cosmic microwave background radiation (CMB), baryon acoustic oscillations (BAG), and supernovae la (SNe-Ia) 

support a constant dark energy equation of state with wq -1. Measuring the evolution of w along the redshift is one of the most 

demanding challenges for observational cosmology. 

Aims. We discuss the existence of a close relation for GRBs, named Combo-relation, based on characteristic parameters of GRB 
phenomenology such as the prompt intrinsic peak energy E^j, the X-ray afterglow, the initial luminosity of the shallow phase Eg, the 
rest-frame duration t of the shallow phase, and the index of the late power-law decay ax. We use it to measure Q.„, and the evolution 
of the dark energy equation of state. We also propose a new calibration method for the same relation, which reduces the dependence 
on SNe la systematics. 

Methods. We have selected a sample of GRBs with 1) a measured redshift z; 2) a determined intrinsic prompt peak energy E^j, and 3) a 
good coverage of the observed (0.310) keV afterglow light curves. The fitting technique of the rest-frame (0.310) keV luminosity light 
curves represents the core of the Combo-relation. We separate the early steep decay, considered a part of the prompt emission, from 
the X-ray afterglow additional component. Data with the largest positive residual, identified as flares, are automatically eliminated 
until the p-value of the fit becomes greater than 0.3. 

Results. We strongly minimize the dependency of the Combo-GRB calibration on SNe la. We also measure a small extra-Poissonian 
scatter of the Combo-relation, which allows us to infer from GRBs alone Q.m = 0.29)^q j5 (Icr) for the ACDM cosmological model, 
and Dm = , 5, Wq = -1.43);|5 gg for the Hat-Universe variable equation of state case. 

Conclusions. In view of the increasing size of the GRB database, thanks to future missions, the Combo-relation is a promising tool for 
measuring D„ with an accuracy comparable to that exhibited by SNe la, and to investigate the dark energy contribution and evolution 
up to z ~ 10. 

Key words. Cosmology: observations. Gamma-ray burst: general. Cosmology: dark energy 


1. Introduction 

Gamma-ray bursts (GRBs) are observed in a wide range of spec¬ 
troscopic and photometric redshifts, up to z ~ 9 (Salvaterra et al. 
2009; Tanvir et al. 2009; Cucchiara et al. 2011). This suggests 
that GRBs can be used to probe the high-z universe, in terms 
of investigating the re-ionization era, population III stars, the 
metallicity of the circumburst medium, the faint-end of galax¬ 
ies luminosity evolution (D’Elia et al. 2007; Robertson & Ellis 
2012; Macpherson et al. 2013; Trent! et al. 2013, 2015), and us¬ 
ing GRBs as cosmological rulers (e.g. Ghirlanda et al. 2004; Dai 
et al. 2004; Amati & Della Valle 2013). The latter works show 
that GRBs, through the correlation between radiated energy or 
luminosity and the photon energy at which their vF(v) spec¬ 
trum peaks Epj, admittedly with a lower level of accuracy, pro¬ 
vide results in agreement with supemovae la (SN la, Perlmutter 
et al. 1998, 1999; Schmidt et al. 1998; Riess et al. 1998), bary- 
onic acoustic oscillations (BAG, Blake et al. 2011), and cosmic 
microwave background (CMB) radiation (Planck Collaboration 
et al. 2013). The universe is spatially flat (e.g. de Bernardis et al. 
2000), and it is dominated by a still unknown vacuum energy, 
usually called dark energy, which is responsible for the observed 
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acceleration. Measuring the equation of state (EOS, a> - pip, 
with p the pressure and p the density of the dark energy) is 
one of the most difficult tasks in observational cosmology to¬ 
day. Current data (Suzuki et al. 2012; Planck Collaboration et al. 
2015) suggest that wo ~ -1 and Wa ~ 0, the expected values 
for the cosmological constant. Although these results are prob¬ 
ably sufficient to exclude a very rapid evolution of dark energy 
with z, we cannot yet exclude that it may evolve with time, as 
originally proposed by Bronstein (1933). In principle, with SNe 
la we can push our investigation up to z ~ 1.7 (Suzuki et al. 

2012) . However we note that the possibility of detecting SNe 
la at higher redshifts will depend on the availability of next- 
generation telescopes and also on the time delay distribution of 
SNe la (see Eigure 8 in Mannucci et al. 2006). 

A solution to this problem may be provided by GRBs: their 
redshift distribution peaks around z ~ 2 - 2.5 (Coward et al. 

2013) and extends up to the photometric redshift of z = 9.4, 
(Cucchiara et al. 2011). Therefore, given this broad range of 
z and their very high luminosities, GRBs are a class of ob¬ 
jects suitable to exploring the trend of dark energy density with 
time (Lloyd & Petrosian 1999; Ramirez-Ruiz & Eenimore 2000; 
Reichart et al. 2001; Norris et al. 2000; Amati et al. 2002; 
Ghirlanda et al. 2004; Dai et al. 2004; Yonetoku et al. 2004; 
Eirmani et al. 2006; Liang & Zhang 2006; Schaefer 2007; Amati 
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et al. 2008; Capozziello & Izzo 2008; Dainotti et al. 2008; 
Tsutsui et al. 2009; Wei et al. 2014). There are two complications 
connected with these approaches. First, the correlations are al¬ 
ways calibrated by using the entire range of SNe la up to z = 1.7; 
therefore, this procedure strongly biases the GRB cosmology 
and for our purposes we need the highest level of independent 
calibration possible and second, the data scatter of these corre¬ 
lations is not tight enough to constrain cosmological parameters 
(in this work we refer to the different extra-scatters published in 
Margutti et al. 2013), even when a large GRB dataset is used. In 
this work we present a method that can override the latter and 
minimize the former. 

Very recently Bernardini et al. (2012) and Margutti et al. 
(2013) (hereafter B12 and M13, respectively) have published an 
interesting correlation that connects the prompt and the after¬ 
glow emission of GRBs (see Fig. 1). This relation strongly links 
the X-ray and y-ray isotropic energy with the intrinsic peak en¬ 
ergy, Ex,iso unveiling an interesting con¬ 

nection between the early, more energetic component of GRBs 
{Ey iso and Ep i) and their late emission (Ex^iso), opportunely fil¬ 
tered for flaring activity. In addition, the intrinsic extra scatter 
is very small cte^jso = 0.30 + 0.06 (Icr). Starting from the re¬ 
sults obtained by B12 and Ml3, and combining them with the 
well-known Epj - Ei^o relation (Amati et al. 2002), we present 
here a new correlation involving four GRB parameters, which 
we named the Combo-relation. This new relation is then used to 
measure Q.,„ and to explore the possible evolution of w of EOS 
with the redshift by using as a “candle” the initial luminosity Lq 
of the shallow phase of the afterglow, which we find to be strictly 
correlated with quantities directly inferred from observations. 

The paper is organized as follows. In Sect. 2 we describe 
in detail how we obtain the formulation of the newly proposed 
GRB correlation, and present the sample of GRBs used to test 
our results, the procedure for the fitting of the X-ray afterglow 
light curves, and finally, the existence of the correlation assum¬ 
ing a standard cosmological scenario. Then, in Sec. 3, we dis¬ 
cuss the use of the relation as cosmological parameter estimator, 
which involves a calibration technique that does not require the 
use of the entire sample of SNe la, as was often done in previous 
similar works. In Sec. 4, we test the use of GRBs estimating the 
main cosmological parameters, as well as the evolution of the 
dark energy EOS. Einally, we discuss the final results in the last 
section. 


2. The Combo-relation 

We present here a new GRB correlation, the Combo-relation, 
obtained after combining Eyjso-Ex,iso~Epj (B12 and Ml3), the 
Eyjso~Epi (Amati et al. 2002) correlations, and the analytical 
formulation of the X-ray afterglow component given in Ruffini 
et al. (2014) (hereafter R14). 

The three-parameter scaling law reported in B12 and M13 
can be generally written as 


where Ex./so is the isotropic energy of a GRB afterglow in the 
rest-frame (0.3-30) keV energy range obtained by integrating 
the light curve in luminosity over a specified time interval, Eyj^o 
is the isotropic energy of a GRB prompt emission, and Epj is the 
intrinsic spectral peak energy of a GRB. Since Ex,iso and £y,,jo 
are both cosmological-dependent quantities, we reformulate this 



Log(E.),,iso/erg)-0,60Log(Epk/l<eV) 

Fig. 1. The Exjso-Eyjso-Ep.i correlation proposed by Bernardini et al. 
(2012) and Margutti et al. (2013) (courtesy R. Margutti). 

correlation in order to involve only one cosmological-dependent 
quantity. 

The right term of Eq. 1 can be rewritten using the well- 
known formulation of the Amati relation, Eyjso ^ which 
provides 

Ex,isc = A EIp (2) 

where j - rj x (5 - 6 and A is the normalization constant. Since 
the Amati relation is valid only for long bursts, in the follow¬ 
ing discussion we will exclude short GRBs with a rest-frame 
Ego duration smaller than 2 s, and short bursts with “extended 
emission” (Norris & Bonnell 2006), or “disguised short” GRBs 
(Bernardini et al. 2007; Caito et al. 2010), which have hybrid 
characteristic between short and long bursts. 

To rewrite Eq. 2 for cosmological purposes, we have calcu¬ 
lated the total isotropic X-ray energy in the rest-frame (0.3-10) 
keV energy range by integrating the X-ray luminosity E(f'), ex¬ 
pressed as a function of the cosmological rest-frame arrival time 
f^, over the time interval (f^ This luminosity E(f^) is ob¬ 

tained by considering four steps (see e.g. Appendix A and Pisani 
etal. 2013). 

It is well known that the X-ray afterglow phenomenology 
can be described by the presence of an additional component 
emerging from the soft X-ray steep decay of the GRB prompt 
emission, and characterized by a first shallow emission, usu¬ 
ally named the plateau, and a late power-law decay behaviour 
(Nousek et al. 2006; Zhang et al. 2006; Willingale et al. 2007). In 
addition, many GRB X-ray light curves are characterized by the 
presence of large, late-time flares, whose origin is very likely as¬ 
sociated with late-time activity of the internal engine (Margutti 
et al. 2010). Since their luminosities are much lower than the 
prompt one, we exclude X-ray flares from our analysis via a 
light-curve fitting algorithm, which will be explained later in the 
text. We then make the assumption that the Ex,,so quantity refers 
only to the component whose X-ray luminosity E(f') is given by 
the phenomenological function defined in R14 

( t' 

E(0 = Eo(i + ^J , (3) 

where Lq, t, and ax are, respectively, the luminosity at f' = 0, 
the characteristic timescale of the end of the shallow phase, and 
the late power-law decay index of a GRB afterglow. 
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Therefore, if we extend the integration time interval to f' j ^ 
0 and f' 2 ^ integral of the function L{t'J in Eq. 3 gives 

Ex, iso - I E{t') dt' = Lq 7t—-[ , (4) 

Jo ll+axi 

with the requirement that ax < -1. This condition is necessary 
to exclude divergent values of Ex,iso computed from Eq. (4), for 
t '^2 worth noting that light curves providing values 

ax > -I could have a change in slope at very late times (beyond 
the XRT time coverage) and/or, in principle, could be polluted 
by a late flaring activity, resulting in a less steep late decay. 

Considering Eqs. 2 and 4, we can Anally formulate the fol¬ 
lowing relation between GRB observables 


La = AEl 


^■'\|l+axl 


-1 


(5) 


which we name the Combo-relation. At first sight, Eq. 5 suggests 
the existence of a physical connection between specific physical 
properties of the afterglow and prompt emission in long GRBs. 
A two-dimensional fashion of the correlation in Eq. 5 in loga¬ 
rithmic units can be written as 


log 



= log 



+ r 


log 



-log 

y 


t/s y 

\l+ax\j 


(6) 

where the set of parameters in the square brackets has the mean¬ 
ing of a logarithmic independent coordinate, and in the following 
will be expressed as log ^X{y, Ep^i, t, ax)]- 

We have tested the reliability of Eqs. 5 and 6 by building a 
sample of GRBs satisfying the following restrictions; 


- a measured redshift z; 

- a determined prompt emission spectral peak energy Epj-, 

- a complete monitoring of the GRB X-ray afterglow light 
curve from the early decay (f' ~ 100 s, when present) un¬ 
til late emission (f^ ~ 10^-10® s). 


We start the analysis by computing the rest-frame (0.3-10) 
keV energy band light curves (see Appendix A). The fitting of 
the continuum part of the X-ray light curves was performed 
by using a semi-automated procedure based on the statis¬ 
tic, which eliminates the flaring part (see e.g. Appendix B and 
Margutti et al. 2011; Zaninoni 2013, and M13, for details). A 
total of 60 GRBs are found, whose distribution in the Combo- 
relation plane is shown in Eig. 2, where the value of the luminos¬ 
ity Lq for each GRB is calculated from the flat ACDM scenario. 
The corresponding best-fit parameters, as well as the extra¬ 
scatter term (Text, have been derived by following the procedure 
by D’Agostini (2005), and are respectively log [A/(erg/s)] = 
49.94 ± 0.27, r = 0.74 ± 0.10, and cr„, =0.33+0.04. The 
Spearman’s rank correlation coefficient is pg - 0.92, while the 
p-value computed from the two-sided Student’s f-distribution, is 
= 9.13 X 10-22. 


3. Calibration of the Combo-relation 

The lack of very nearby (z ~ 0.01) GRBs prevents the possibility 
of calibrating GRB correlations, as is usually done with SNe la. 
In recent years different methods have been proposed to avoid 
this “circularity problem” in calibrating GRB relations (Kodama 
et al. 2008; Demianski et al. 2012, and references therein). The 



Log[^^,,/keV] -1.36 X Log[(T/s)/|l+a|] 

Fig. 2. Plot of the correlation considering the entire sample of 60 GRBs. 
The green empty boxes are the data of each of the sources, derived 
as described in Appendix B, the solid black line is the best fit of the 
data, while the dotted grey lines and the dashed grey lines correspond, 
respectively, to the dispersion on the correlation at IcTex and Str^x- 



redshift z 


Fig. 3. The residual distance modulus - p,* for different values of 

the density cosmological parameters and Da up to z = 2.0. We 

consider the best fit to be the standard ACDM model, where = 0.27, 
Da = 0.73, and Hf, = 1\ km/s/Mpc (black line). Union2 SNe la data 
residuals are shown in grey. The large spread (more than 1 magnitude) 
shown by p at z = 1.5 and at z = 0.145 (the two vertical dashed lines) 
where the scatter is almost 0.2 magnitudes is clearly evident. 


common approach uses an interpolating function for the distri¬ 
bution of luminosity distances (distance moduli) of SNe la with 
redshift, and then extending it to GRBs at higher redshifts. In the 
following we introduce and describe an alternative method for 
calibrating the Combo-relation (see also Ghirlanda et al. 2006) 
which consists of two steps: 1) we identify a small but sufficient 
subsample of GRBs that lie at the same redshift, and then infer 
the slope parameter y from a best-fit procedure; 2) once we de¬ 
termine y, the luminosity parameter A can be obtained from a 
direct comparison between the nearest, z = 0.145, GRBs in our 
sample and SNe la located at very similar redshifts. We note that 
this approach is different from previous ones in that we do not 
use the whole redshift range covered by SNe la, but we limit our 
calibration analysis to z = 0.145 where the effect of the cosmol¬ 
ogy on the distance modulus of the calibrating SNe la is small 
(see Fig. 3). 
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Table 1. List of the five GRBs used for the determination of the slope parameter 7 of the correlation. In the Col. 1 is shown the GRB name, in Col. 
2 the redshift z, in Col. 3 the intrinsic peak energy of the burst, in Col. 4 the initial flux of the afterglow additional component Fg in the common 
rest-frame energy range (0.3 - 10) keV, in Col. 5 the time parameter (t), in Col. 6 the late power-law decay index (ax)- 


GRB 

z 

Ep, (keV) 

Fg (erg cm ') 

t(s) 

ax 

060729 

0.54 

104 ± 24 

(2.32 ± 0.05) X 10-“ 

61462 ± 2844 

-1.71 ±0.03 

081007 

0.5295 

61 ± 15 

(4.48 ± 0.64) X lO^" 

2178 ±457 

-1.17 ±0.04 

090424 

0.544 

273 ±5 

(4.20 ± 0.26) X 10-'* 

204 ± 16 

- 1.20 ± 0.01 

090618 

0.54 

257 ± 41 

(1.18 ± 0.03) X 10-'* 

1206 ±47 

-1.46 ±0.01 

100621A 

0.542 

146 ± 23 

(1.75 ± 0.22) X 10^'° 

3613 ± 746 

-1.05 ±0.03 



Fig. 4. Plot of the correlation found for the sample of five GRBs lo¬ 
cated at the same redshift. The blue triangles are the data of each of 
the five sources, derived from the procedure in Appendix B. The solid 
line represents the best fit while the dashed line is the dispersion on the 
correlation at Icr^r- 


3.1. The determination of the siope y 

The existence of a subsample with a sufficient number of GRBs 
lying at almost the same cosmological distance would, in prin¬ 
ciple, allow us to infer 7 , overriding any possible cosmological 
dependence (assuming a homogeneous and isotropic universe). 
In our sample of 60 GRBs there is a small subsample of 5 GRBs 
located at a very similar redshift, see Table 1. The difference be¬ 
tween the maximum redshift of this 5-GRB sample (z - 0.544) 
and the minimum one (z - 0.5295) corresponds to a variation 
of 0.015 in redshift and of 0.07 in distance modulus p in the 
case of the standard ACDM model. This very small difference is 
sufficient for our purposes. 

To avoid any possible cosmological contamination, we do 
not consider the luminosity Lq as the dependent variable, but we 
consider instead the energy flux Fq, which is related to the lu¬ 
minosity through the expression Lq = ^o- This assumption 

does not influence the final result, since di is almost the same for 
the 5-GRB sample and, therefore, the Andf term can be absorbed 
in the normalization constant. Consequently, we build the energy 
flux light curve for each GRB, and then, following the same pro¬ 
cedure described in Sec. 2, we perform a best-fit analysis of this 
subsample of five GRBs using the maximum likelihood tech¬ 
nique. From the best fit we find a value of 7 = 0.89 + 0.15, with 
an extra scatter of (Text - 0.40 + 0.04, see Fig. 4. 

3.2. Caiibration of A with SNe la 


Table 2. List of the 5 SNe la selected from the Union2 sample 
(Amanullah et al. 2010; Suzuki et al. 2012) and used for the calibration 
of the GRB correlation. In the first column it is shown the name of the 
SN, in the second column the redshift and in the third the distance mod¬ 
ulus obtained from the light curve fitting with the SALT2 SED method 
(Guy et al. 2007). 


SN 

z 

F 

1999bm 

0.1441 

38.836 ±0.157 

10106 

0.14629 

39.559 ±0.120 

20051n 

0.14567 

39.057 ±0.126 

2005gx 

0.144621 

39.291 ±0.113 

20051d 

0.1437 

39.186 ±0.116 


is clear that at these redshifts (z — 0.145) the offset provided by 
SNe la is much smaller than the value inferred from SNe la at 
larger redshifts, see Fig. 3. In this light we have selected Union2 
SNe la (Amanullah et al. 2010; Suzuki et al. 2012) with red¬ 
shift between z = 0.143 and z = 0.147. We find five SNe la that 
satisfy this condition; their properties are shown in Table 2. We 
then compute the average value of their distance modulus and 
its uncertainty, (p) = 39.19 + 0.27, which will be used for the 
final calibration. The parameter A can be obtained inverting Eq. 
6 , and considering Lq = Andfi Fq: 




-H log 47r -I- log 


Fq \ ^ 
erg/cm^/s / 


- ylogXiy,Ep,i,T,ax). 


(7) 


The generic SN distance modulus p can be directly related to the 
luminosity distance di by 


p = 25 + 5 log 


di 

Mpc 


= -97.45-H 5 log I-^'1 , 
cm / 


( 8 ) 


where the last equality takes into account the fact that di is ex¬ 
pressed in cm. Substituting this last expression for log (di) in Eq. 
7 we obtain 
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(ip) + 97.45) + if/(y, Epj, t, ax, Fq) , 


(9) 


where the term Ep i, t, ax, Fq) comprises the last three terms 
on the right hand side of Eq. 7. Substituting the quantities for 
GRB 130702A in Eq. 9, and the value of (p) previously obtained, 
we infer a value for the parameter log[A/(erg/s)] = 49.54 + 0.20, 
where the uncertainty also takes the cr„, value found above into 
account. 


4. Cosmology with the Combo-relation 

4. 1. Building the GRB Hubble diagram 


Among the considered 60 GRBs, the nearest one is GRB 
130702A at z = 0.145. We use it to perform the calibration of 
the Combo-relation with SNe la located at the same distance. It 


We now discuss the possible use of the proposed GRB Combo- 
relation to measure the cosmological constant and the mass den¬ 
sity, as well as their evolution with redshift z. The possibility of 
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estimating the luminosity distance di from the GRB observable 
quantities allows to us dehne a distance modulus for GRBs, and 
then its uncertainty, as 

^tCRB = -97.45 + ^ [54 - Epj, t, ax, T^o)], (10) 

where = log[A/(erg/s)]. The quantity ficRs can be directly 
compared with the theoretical cosmological expected value p,/,, 
which depends on the density parameters and the curva¬ 
ture term Ga = 1 - - Ga, and the Hubble constant Hq 


Eih = 25 -H 5 log di(z, Q,„, Ga, Hq) . 
The luminosity distance di is dehned as 

-J— sinh r _ 

Vind [Jo 

r 

Jo 


d,^—(l+z) 

Hq 


V£(z'.n„.nA,wfe)) 


dz' 


y/E{z',0.m,O.A,wiz)) 


1 

I vra 


sm 


r 


VindA-' 




( 11 ) 

Ga >0 

( 12 ) 

Ga <0 


with E{z, nm, Q.A, w{z)) = G„,(l -hz)3 -hGa(1 -hQa(1 + 

z)^ (see e.g. Goobar & Perlmutter 1995). In the following we 
fix the Hubble constant at the recent value inferred from low- 
redshift SNe la, corrected for star formation bias, and calibrated 
with the LMC distance (Rigault et al. 2014); Hq = 70.6 + 2.6. 

The corresponding uncertainties on /ucrb were computed 
considering an observed term, cr^obs, which takes into account 
each uncertainty on the observed quantities of the Combo- 
relation, e.g. Fq, t, a, and Epj, and a “statistical” term, crfirei, 
which takes into account the uncertainties on the parameters of 
the Combo-relation, A and y, and the weight of the extra scatter 
value (Text- The hnal uncertainty on each single GRB distance 
modulus 

O-flCRB = ^{o-Hobs9 + (cTllrel)^ ( 13 ) 

allows us to build the Combo-GRB Hubble diagram (see also 
Izzo et al. 2009), which is shown in Fig. 5. It is possible to quan¬ 
tify the reliability of any cosmological model with our sample 
of 60 GRBs, which represents a unique dataset from low red- 
shift (z = 0.145) to very large distances (z = 8.23). To reach our 
goal, we make the fundamental assumption that our GRB sam¬ 
ple is normally distributed around the best-fit cosmology, which 
we are going to estimate. With this hypothesis, we consider as 
test-statistic the chi-square test, which is defined as 




^ [EcrbAz, - Ethiz, Qa, w(z))f 


^EcRB.i 


(14) 


where hcrb{z, Jl, (p), ^ithiz, G„,, Ga, vv(z)), and ct^icrb are respec¬ 
tively dehned in Eqs. 10, 11, and 13. The cosmology is included 
in the quantity ;U,/,(z, G,„, Hq), which we allow to vary. To 
determine the best conhguration of parameters that most closely 
hts the distribution of GRBs in the Hubble diagram we maxi¬ 
mize the log-likelihood function, -2 ln(e^ ), which is equivalent 
to the minimization of the function dehned in Eq. 14. 


4.2. Fit results 
4.2.1. ACDM case 



Fig. 5. The Combo-GRB Hubble diagram. The black line represents the 
best ht for the function /j(z) as obtained by using only GRBs and for the 
case of the ACDM scenario. 



Fig. 6. 1-cr (Ay^ = 2.3) conhdence region in the (Hn,, Da) plane (left 
panel) and in the (On,, Wq) plane (right panel) for the Combo-GRB sam¬ 
ple (dark blue), and for the total (60 observed -l- 300 MC simulated 
GRBs, (light blue). 


4.2.2. Variable wq case 


In a Hat Universe (G = -H Ga = 1) with a constant value 
of wq different from the standard value w - -1, we can pro¬ 
vide useful constraints for alternative dark energy theories. In 
this case, we only vary the matter density and the dark energy 
equations of state, obtaining an estimate of the density mat¬ 
ter of Q„, = 0.40;(Qjg and of a dark energy EOS parameter 


Wo = -1.52 


+ 0.94 

- 0 . 93 ' 


In the ACDMmodel, the energy function E(z,Ei^,Q.A,w{z)) is 
characterized by an EOS for the dark energy term hxed to w = 
-1. Since we have that 0,„ ■¥ Qa -h Qa = 1, we vary the matter 
and cosmological constant density parameters, also obtaining in 
this way an estimate of the curvature term. We obtain that GRBs 
alone provide Q,„ = 0.29 ^qj 5, see also Eig. 6. 


4.2.3. Evolution of w(z) 

An interesting case-study consists of a time-evolving dark en¬ 
ergy EOS in a hat cosmology, since the evolution of w(z) can 
be directly studied with GRBs at larger redshifts. We consider 
an analytical formulation for the evolution of w with the red- 
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Fig. 7. l-cr - 2-3) confidence region in the (fin,, Wg) plane (right 
panel) for the Combo-GRB sample (dark blue), and for the total of 60 
observed + 300 MC simulated GRBs (light blue). The black dashed line 
represents the l-cr confidence region obtained using the recent Union 
2.1 SNe la sample (Suzuki et al. 2012). 


where ^ - z, Epj, t, and |aH-1|, and we found the following mean 
values and dispersions: /U^+ctj = 0.26+0.27,= 2.54+ 
0.40, Pt±o-t - 2.85 + 0.65, and p\a+\\ + = -0.36 + 0.31. 

Then, from these distributions we computed \ogX(y,Ep,i,T,a) 
and we generated the initial luminosity log Lq = T log X + A 
from the frequency distribution 


/(logic) = -exp 

ylnCText 


(logip -ylogX-A)^ 


(17) 


assuming that the Combo-relation is independent of the redshift 
and considering its extra-scatter cr„, . The values of y. A, and 
(Text are reported in Section 3. Finally, to complete the set of pa¬ 
rameters necessary to compute the distance modulus of the sim¬ 
ulated sample of GRBs from each pair (logic, z), we generated 
the corresponding log Fq (p.p^ + (T/r„ = -9.87 + 0.85). In the fol¬ 
lowing, the attached errors on the MC simulated Combo-relation 
parameters will be taken as 30% of their corresponding values, 
which reflects the uncertainty of the “real” GRB sample. 

We have verified whether the constraints on (G,^, Q^) will 
improve by using a larger sample of 360 GRBs (the real sam¬ 
ple of 60 GRBs observed and a MC-simulated sample of 300 
GRBs, described above). The improvement on the constraints 
on Qtn and Qa is clear: the uncertainties on the density parame¬ 
ters improve considerably (Q„, = 0.27^® °® for the ACDM case, 
wc = - 1 . 16 );Q 3 g for the variable wq case), as is also clear in the 
contour plots shown in Figs. 6 and 7. 


shift, which was proposed by Chevallier & Polarski (2001) and 
Linder (2003) (CPL), and where the w{z) can be parameterized 
by 


z 

w(z)^Wo + Wa- -. (15) 

1 +z 

The CPL parameterization implies that for large z the w(z) 
term tends to the asymptotic value wo + Wa- Using a sample ex¬ 
tending at large redshifts, e.g. GRBs, will allow a better estimate 
of these parameters since the effects of a varying w(z) on the dis¬ 
tance modulus are more evident at redshift z > 1. The GRB data 
provide a best-fit result (wc = -1.43;!;° ^^, = L87;|;2 57 )' 

4.2.4. Perspective: a Monte Carlo simulated sample of 
GRBs 

The current sample of GRBs that satisfies the Combo-relation is 
quite limited (60 bursts), when compared to the sample in Amati 
& Della Valle (2013) (~ 200 bursts) or the SNe sample in the 
Union 2.1 release (Suzuki et al. 2012). A more numerous sample 
can help to understand whether the Combo-relation can provide 
better constraints on the cosmological parameters. To this aim, 
following the prescription of Li (2007), we used Monte Carlo 
(MC) simulations to generate a sample of 300 synthetic GRBs 
satisfying the Combo-relation. This value comes from the ex¬ 
pected number of GRBs detected in five years of operations of 
current (Swift) and future (SVOM (Gotz et al. 2009), and LOFT 
(Feroci et al. 2012)) missions dedicated to observing GRBs. 

First, we fitted the log-normal distributions of the 60 ob¬ 
served z, Ep i, T, and \a + 1 |, 


/(log^) = 




exp 




2cr| 


(16) 


4.3. GRBs compared with SNe, CMBs, BAOs 

In order to compare the Combo-GRB sample results, we also 
consider the following datasets: 

- the measurements of the baryon acoustic peaks Ao^^ = 
(0.474 + 0.034,0.442 + 0.020,0.424 + 0.021) at the corre¬ 
sponding redshifts zbao - (0.44,0.6,0.73) in the galaxy cor¬ 
relation function as obtained by the WiggleZ dark energy 
Survey (Blake et al. 2011). The BAG peak is defined as 
(Eisenstein et al. 2005) 


^BAO - 


czbao rizBAO, Elm, El A, w(z))^ 


Hq E{z, Elm, EIa, w(z)) 


1/3 






CZbAO 


(18) 


where r(z. Elm, EIa, w(z)) is the comoving distance. The best- 
fit cosmological model is determined by the minimization of 
the corresponding chi-squared quantity 


3 

xIbAO - '^i^BAO - A(obs,i)YC ^{AbaO - A(obs,i)), (19) 

i=l 


where C ' is the inverse covariance matrix of the measure¬ 
ments of the WiggleZ survey (Blake et al. 2011). 

- the measurement of the shift parameter Robs - 1.7407 + 
0.0094 as obtained from the Planck first data release (Planck 
Collaboration et al. 2013). The R quantity is the least cosmo¬ 
logical model-dependent parameter (particularly from Hq) 
that can be extracted from the analysis of the CMB (Wang 
& Mukherjee 2006) and is defined as 


R = 



dz' 

E(z',Elm,ElA,w(z))’ 


( 20 ) 
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Fig. 8. 1-cr = 2.3) confidence region in the Da) plane for the 
observed GRB sample (blue), with the inclusion of the MC-simulated 
300 GRBs sample (cyan), and with the samples of SNe (grey), BAOs 
(red), and CMBs (green). The dashed line represents the condition of 
the Flat Universe flm + Qa = 1- 


where Zrec is the redshift of the recombination. The best-fit 
cosmological model is determined by the minimization of 
the corresponding chi-squared term 


XcMB ~ 


(Robs - Rf 


crR 


2 

obs 


( 21 ) 


We use a grid-search technique to vary the values of the cos¬ 
mological which are used to solve numerically the chi-squared 
equations defined above. Every dataset provides the respective 
distribution of the cosmological parameters, so when combining 
datasets we simply determine the best-fit model as the sum of the 
single x^i values for any pair of parameters and for the combined 
values, 

2 2 2 2 2 
Xtot - XcRB XsNe XbAO XcMB^ ( 22 ) 

and make the final plots using the Mathematical software suite. 
Cosmological parameter uncertainties were estimated from sin¬ 
gle and combinedstatistics for each dataset. In Fig. 8 we show 
the parameter spaces with the contour at the l-cr (Ax^ - 2.3) 
confidence limit for each considered dataset, and for the ACDM 
case. We note that in all the datasets and fits, the Hubble constant 
is fixed to the value of Hq - 70.6+2.6 (Rigault et al. 2014). From 
the combination of all the considered datasets, for the ACDM 
case, we obtain a matter density value of Q„, = 0.28;(yy[, which 
shows no strain with the Combo-GRB results and the simulated 
one, see Fig. 8. 


5. Discussions and conciusions 

In this paper we have presented the “Combo” relation, a new 
tool for GRB cosmology. This relationship provides a very close 

' http://www.wolfram.com/mathematica/ 


link between prompt and afterglow parameters and it is charac¬ 
terized by a small intrinsic scatter, which makes this correlation 
very suitable for cosmological purposes. We recognize the fun¬ 
damental role of the Swift-XRI (Gehrels et al. 2004; Burrows 
et al. 2005) which, thanks to its ability to slew very rapidly to¬ 
ward the location of a GRB event, provides real-time and de¬ 
tailed data of GRB afterglow light curves, whose evolution is at 
the base of the proposed Combo-relation. From our analysis the 
following results emerge: 

- the proposed two-step calibration of the Combo-GRB rela¬ 
tion greatly minimizes the dependence on SNe la; 

- GRBs data alone provide for the ACDM case Q.„, = 

0.29!°;23, see Fig. 8; 

- a recent paper (Milne et al. 2015) highlights the existence 
of an observational bias (a systematic difference in the ve¬ 
locity of SNe la ejecta, which is reflected in their curves), 
potentially affecting the measurements of cosmological pa¬ 
rameters obtained with SNe la. On the basis of our results 
we conclude that given the current accuracy of GRB mea¬ 
surements we cannot exclude, within the errors, that an effect 
like this is at play; however, this effect should not change the 
conclusions derived from SNe-Ia observations; 

- the launch of advanced and more sensitive detectors, such 
as the incoming SVOM (Gotz et al. 2009) and the proposed 
LOFT (Feroci et al. 2012) missions (and the expected Swift 
operations in the near future), will dramatically increase the 
number of GRBs in the dataset. In five years of operation 
of the SVOM mission alone, we expect to reach a “good 
enough” sample of 300 GRBs. With a Monte Carlo simu¬ 
lated sample of 300 GRBs, we will significantly improve the 
accuracy of measurement, up to Q„, = 0.27 /qq 5, which 
is comparable with type la SNe (Suzuki et al. 2012); 

- by using the CPL analytical parameterization, adopted to 
study the evolution of the dark energy EOS (see Eq. 15), we 
find = 0.40!0;22, wq = -1-43:°;™, and = 1.87!^;38; 

- the analysis of a combined (SNe+BAO+CMB) dataset con¬ 
firms that the increasing size of the GRB sample will im¬ 
prove the accuracy of the measurement of the Q.„, parameter 
and in particular of the evolution of w up to z ~ 10. 

- the analytical expression of the Combo-relation provides an 
explicit close link, characterized by a small intrinsic scatter, 
between the prompt and the afterglow GRB emissions. This 
points out the existence of a physical connection between the 
prompt and the afterglow emissions, which represents a new 
challenge for GRB models. 
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Appendix A: Computation of the rest-frame 0.3-10 
keV luminosity L{Q 

The rest-frame 0.3-10 keV luminosity L{t'^ was obtained by 
considering four steps. 


(1) We obtained the Swift-XKT flux light curves in the observed 
0.3-10 keV energy band^. 

(2) We transformed the observed flux fgh^ from the observed en¬ 
ergy band 0.3-10 keV to the rest-frame energy band 0.3- 
10 keV by assuming an absorbed power-law function as the 
best bt for the spectral energy distribution of the XRT data, 
N(E) ~ E^'^, with Galactic and intrinsic column densities 
obtained from the H1 radio map (Kalberla et al. 2005) and 
from the best fit of the total afterglow spectrum, respectively. 
By using the photon indexes inferred for each time interval, 
the rest-frame flux light curve frf is given by 


(3) We transformed the observed time G into the rest-frame time 
by correcting for z 

4 = W(l+z)- (A.2) 

(4) We defined the isotropic luminosity as 

L = And^frf. (A.3) 


Appendix B: Determination of the sample and 
verification of the correlation 


To obtain the parameters involved in the Eq. 5, we needed to se¬ 
lect an adequate sample of GRBs, to fit their X-ray light curves, 
and to collect or to calculate their Epj values. The selection cri¬ 
teria have already been delineated in Sec. 2. 

The entire procedure works in the rest-frame 0.3-10 keV en¬ 
ergy range for all GRBs. For the X-ray light curve fitting tech¬ 
nique we developed a semi-automated code performing all the 
needed operations. The code is based on the IDL^ language, and 
the fitting routine used is MPFIT'^ (Markwardt 2009), which is 
based on the non-linear least squares fitting. First, the procedure 
fits the complete light curve, then it eliminates at every iteration 
the data point with the largest positive residual, until it obtains a 
fit with a p-value greater than 0.3. To fit the light curves consid¬ 
ered in luminosity units (erg/s), we use the composite function 
(R14): 

1 . a power law for the early steep decay 


Lit) = Lp 



(B.l) 


with Lp the normalization factor and Up the slope; 
2. Eq. 3 for the afterglow additional component. 


An application of this joint fitting procedure is shown in Fig. 
?? for the case of GRB 060418. After the fitting procedure, we 
select only the GRBs with ax < -1, a condition necessary for 
the convergence of the integral in Eq. 4. The final total sample, 
summarized in Table B.l, consists of 60 GRBs. 


^ http://www.swift.ac.uk/burst_analyser/ 

^ Interactive Data Language, http://www.exelisvis.com/language/en- 
US/ProductsServices/TDL.aspx 
^ http://purl.com/net/mpfit 
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Table B.l. Long bursts with < ”1 analysed in this work (first column) and their main parameters: GRB name (first column), redshift z (second column), intrinsic 
peak energy Ep^i (third column), flux in the 0.3-10 keV rest-frame energy band Fq (fourth column), the time parameter t (fifth column), and the late power-law 
decay index ax (sixth column). All errors are at the l-cr confidence level. 


GRB 

z 

Epj (keV) 

Fo (erg cm ^s ^) 

t(s) 

ax 

050401 

2.9 

467 ± 110 

(2.35 ± 0.13) X 10-"’ 

733 ± 102 

-1.64 ±0.07 

050922C 

2.198 

415± 111 

(3.65 ± 0.39) X lO-’" 

234 ± 35 

-1.51 ±0.04 

051109A 

2.346 

539 ± 200 

(2.24 ± 0.29) X lO-’" 

429 ± 62 

-1.26 ±0.02 

060115 

3.53 

285 ± 34 

(5.07 ± 0.84) X 10-’^ 

12553 ± 6079 

-2.87 ±0.72 

060418 

1.489 

572 ±143 

(6.00 ± 1.10) X lO-’” 

294 ± 62 

-1.65 ±0.06 

060526 

3.21 

105 ±21 

(5.05 ± 0.49) X 10- 

25847 ± 9019 

-4.44 ± 0.97 

060707 

3.43 

279 ± 28 

(1.94 ± 0.53) X 10-" 

1067 ± 485 

-1.11 ±0.06 

060729 

0.54 

104 ± 24 

(2.28 ± 0.05) X 10-" 

61462 ± 2844 

-1.71 ±0.03 

060814 

0.84 

473 ± 155 

(5.04 ± 0.75) X 10-” 

3183 ±682 

-1.30 ±0.05 

060908 

1.8836 

514±102 

(1.02 ± 0.25) X lO-"* 

207 ± 72 

-1.67 ±0.12 

060927 

5.467 

475 ± 47 

(1.74 ± 0.23) X 10-" 

230 ± 73 

-1.44 ±0.16 

061121 

1.314 

1289 ± 153 

(3.08 ± 0.10) X lO-’" 

1544 ± 83 

-1.54 ±0.02 

071020 

2.145 

1013 ± 160 

(3.50 ± 1.30) X lO-"* 

9.9 ± 4.0 

-1.16 ±0.02 

080319B 

0.937 

1261± 65 

(6.52 ± 0.37) X 10-'’ 

22.11 ±0.97 

-1.70 ±0.01 

080413B 

1.1 

150 ±30 

(4.01 ±0.35) X lO-’" 

184 ±24 

-1.11 ±0.02 

080605 

1.6398 

650 ± 55 

(3.48 ± 0.16) X lO-"* 

83.7 ± 6.0 

-1.42 ±0.02 

080607 

3.036 

1691 ±226 

(2.31 ±0.54) X lO-’" 

387 ± 91 

-1.65 ±0.06 

080721 

2.591 

1741± 227 

(2.61 ±0.13) X lO-"* 

431±100 

-1.61 ±0.02 

080810 

3.35 

1470 ± 180 

(1.27 ± 0.50) X lO-’" 

647 ± 239 

-1.83 ±0.10 

080916A 

0.689 

184 ± 18 

(5.05 ± 0.84) X 10-" 

967 ±251 

-1.06 ±0.03 

080928 

1.692 

95 ±23 

(3.83 ± 0.56) X 10-" 

5976 ± 1499 

-3.11 ±0.33 

081007 

0.5295 

61 ± 15 

(4.61 ±0.66) X 10-" 

2178 ±457 

-1.17 ±0.04 

081008 

1.9685 

261 ± 52 

(6.46 ± 0.79) X 10-" 

1666 ± 368 

-1.86±0.13 

081222 

2.77 

505 ± 34 

(1.58 ± 0.18) X lO-’" 

699 ± 93 

-1.61 ±0.04 

090102 

1.547 

1149± 166 

(5.00 ± 1.40) X lO-"* 

66 ± 18 

-1.43 ±0.02 

090418A 

1.608 

1567 ± 384 

(4.23 ± 0.26) X lO-’" 

553 ± 64 

-1.63 ±0.06 

090423 

8.2 

437 ± 55 

(6.01 ± 0.75) X 10-’2 

376 ± 115 

-1.19±0.11 

090424 

0.544 

273 ±5 

(4.32 ± 0.27) X lO-"* 

204 ± 16 

-1.20 ±0.01 

090516 

4.109 

971 ± 390 

(4.10 ± 1.40) X 10-" 

761 ±382 

-1.38±0.12 

090618 

0.54 

257 ± 41 

(1.22 ± 0.03) X lO-"* 

1206 ± 47 

-1.46 ±0.01 

091018 

0.971 

55 ±20 

(8.25 ± 0.56) X 10-"* 

169 ± 20 

-1.28 ±0.02 

091020 

1.71 

280 ± 190 

(3.94 ± 0.33) X 10-"* 

311 ±36 

-1.36 ±0.03 

091029 

2.752 

230 ± 66 

(1.39 ± 0.08) X 10-" 

2806 ± 343 

-1.31 ±0.04 

091127 

0.49 

53.6 ±3.0 

(1.03 ± 0.07) X lO-** 

1966 ± 165 

-1.51 ±0.02 

100621A 

0.542 

146 ± 23 

(1.81 ±0.22) X 10-"* 

3613 ±746 

-1.05 ±0.03 

100814A 

1.44 

259 ± 34 

(1.72 ± 0.08) X 10-" 

36217 ±4018 

-2.00 ± 0.08 

100906A 

1.727 

289 ± 46 

(1.25 ± 0.17) X 10-"* 

2341 ± 402 

-2.19±0.11 

110213A 

1.46 

231 ±21 

(4.09 ± 0.31) X 10-"* 

2037 ± 209 

-2.13 ±0.07 

110422A 

1.77 

421 ± 14 

(8.40 ± 0.98) X 10-"* 

358 ± 56 

-1.41 ±0.03 

111228A 

0.714 

58.3 ±5.1 

(7.66 ± 0.53) X 10-" 

4094 ± 464 

-1.36 ±0.03 

120119A 

1.728 

515 ±22 

(2.00 ± 0.28) X 10-" 

8597 ± 2226 

-2.57 ±0.26 

120811C 

2.671 

198 ± 11 

(7.20 ± 0.63) X 10-" 

1204 ± 340 

-1.57±0.16 

120907A 

0.97 

303 ± 65 

(7.26 ± 0.88) X 10-" 

481 ± 98 

-1.13 ±0.04 

120922A 

3.1 

303 ± 14 

(1.82 ± 0.67) X lO-’" 

119±56 

-1.10 ±0.04 

121128A 

2.2 

198 ± 15 

(2.61 ±0.59) X lO-’" 

1244 ± 363 

-2.25 ±0.18 

121211A 

1.023 

194 ± 26 

(2.61 ±0.56) X 10-" 

2296 ± 893 

-1.23 ±0.10 

130408A 

3.757 

1294 ± 190 

(1.18 ± 0.24) X lO-’" 

216± 118 

-1.32±0.16 

130420A 

1.297 

129 ±7 

(4.90 ± 0.86) X 10-” 

571 ± 158 

-1.05 ±0.03 

130427A 

0.3399 

1161 ±7 

(1.42 ± 0.07) X 10-” 

996 ± 61 

-1.40 ±0.01 

130505A 

2.27 

2063 ± 101 

(4.54 ± 0.44) X lO-’" 

1667 ± 182 

-1.63 ±0.03 

130610A 

2.092 

912 ± 133 

(2.20 ± 1.60) X lO-’" 

41 ±35 

-1.12 ±0.06 

130612A 

2.006 

186±31 

(1.00 ± 0.46) X 10-” 

602 ± 468 

-1.31 ±0.24 

130701A 

1.155 

192 ±9 

(2.14 ± 0.53) X lO-"* 

50± 15 

-1.28 ±0.03 

130702A 

0.145 

14.9 ± 2.3 

(3.53 ± 0.46) X 10-" 

61448 ± 10744 

-1.35 ±0.03 

130831A 

0.4791 

81 ±6 

(1.87 ± 0.29) X lO-’" 

1925 ± 332 

-1.73 ±0.07 

131030A 

1.295 

406 ± 22 

(8.08 ± 0.49) X lO-’" 

370 ± 41 

-1.31 ±0.02 

131105A 

1.686 

548 ± 83 

(3.96 ± 0.36) X 10-" 

986 ± 235 

-1.19 ±0.07 

131117A 

4.042 

222 ± 37 

(2.04 ± 0.51) X 10-" 

291 ± 115 

-1.32±0.10 

140206A 

2.73 

448 ± 22 

(1.86 ± 0.16) X lO-’” 

1477 ± 154 

-1.52 ±0.03 

140213A 

1.2076 

177 ±4 

(1.84 ± 0.24) X lO-’" 

1870 ±337 

-1.31 ±0.04 


9 



